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Abstract 

Certain physical aspects of quantum error correction are discussed for a 
quantum computer (n-qubit register) in contact with a decohering environ- 
ment. Under rather plausible assumptions upon the form of the computer- 
environment interaction, the efficiency of a general correcting procedure is 
evaluated as a function of the spontaneous-decay duration and the rank of 
errors covered by the procedure. It is proved that the probability of errors 
can be made arbitrarily small by enhancing the correction method, provided 

the decohering interaction is represented by a bounded operator. 
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I. INTRODUCTION 



Since Shor |IJ demonstrated that the classical factoring problem can, in principle, be 
efficiently solved on a quantum computer, a multitude of papers dealing with quantum 
computing have emerged. It was soon recognized BJ3], however, that the advantage of im- 
plementing entangled quantum states may be spoilt by their vulnerability to errors. Namely, 
the destructive interference of the omnipresent environment leads to an exponential loss of 
the probability that the computation runs in the desired way. The number of runs needed 
to perform a successful calculation then increases exponentially which takes one back to the 
purview of complexity problems. 

These pessimistic views have abated to a certain extent since the first work pioneering 
quantum error correction appeared ||. It immediately became a subject of fascination that 
something as delicate to handle as quantum state can be mended without knowing any 
particulars about it. In view of that, a variety of error correcting codes || and related 
complexity requirements have been thoroughly discussed. 

The basic idea of quantum correcting is fairly simple. The computational state is en- 
coded, by introducing redundancy, into a more robust one, which can then be rectified, 
provided only errors from a certain sub-class have occurred. It is subsequently argued that 
the remaining errors appear with a small probability, which is a key prerequisite for the 
proof of correction effectiveness. A typical example of the above sub-class are single-qubit 
errors. In this case, simultaneous errors on more than one qubit are expected to conform 
to the law of independent probabilities, and, therefore, to be dominated by the single-qubit 
errors. If the latter errors can be put away, the system is indeed much less susceptible to 
perturbations. 

While quantum codes and related topics have been well explored during over a year 
of their existence, little has been said on the genuine physical aspects of the quantum 
error correction, though some of them certainly deserve our attention. Namely, despite all 
correction methods resemble the watch-dog stabilization (discussed, e.g., in ||), there is 
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an important difference: For the watch-dog stabilization (in the ideal case) to function, 
corrections have to be repeated at an infinite rate. But, in reality, there is always a principal 
bound on the correction frequency (apart from others, for the quantum computation must 
not be interfered with) — the exponential decay of the original state cannot be avoided. Of 
course, this does not necessarily imply that the error correction brings no profit. In order to 
estimate the correction efficiency, however, the temporal aspects of the correction process 
have to be carefully scrutinized. As shown in the following example, this cannot be done 
irregardless the particulars of the decohering interaction. 

Consider a simple correction scheme capable of eliminating single- and double-qubit 
amplitude errors (i.e., |0) <-> |1) flips — see sect. II): States of a single qubit are encoded into 
states of a qubit pentet, |0) — > |00000) and |1) — > 1 11111) , and the correction is represented 
by the transformation {b^bsb^) —> \bbbbb) on the encoded states, where b is equal to bi that 
occurs at least three times in the initial state (majorization rule). Let us try to compare the 
stabilization effect of the above correction under two different error-producing interactions: 
H^> = Y,i hwi&l an d = J2k^i ^ki^x^xi where a x is the first Pauli matrix (a x : |0) <-> |1)) 
operating on the i-th qubit. As both the interactions cause only amplitude errors, the 
correction procedure is, in principle, applicable. The corresponding evolution operators have 
the form: U^it) = l\i[cos(uJit) +i<J l x sin(c<;^)] and U^ 2 \t) = H kitl [cos(ukit) +ia^a l x sin(ukit)]. 
By breaking the products into parts, it is evident that the error proliferation is much faster 
for the second hamiltonian. Namely, for short times t, the errors which cannot be cured by 
our procedure (i.e., the terms in U(t) corresponding to three- and higher-qubit errors) are of 
order t 3 in the first case and t 2 in the second one. To achieve the same level of stabilization 
for both interactions, the corrections in the second case have to be repeated at much higher 
rate. 

The plan of the paper is as follows. In the second section we briefly recall the basic 
facts on decoherence and quantum error correction. The third section is devoted to the 
evaluation of the correction efficiency, under a certain mild assumption on the form of 
decohering interaction. In the last section we present conditions under which a rigorous 
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meaning can be given to our formulas, and discuss their applicability to realistic situations. 
In particular, we prove that the quantum computer can be stabilized if the decohering 
interaction is represented by a bounded operator. 



II. DECOHERENCE AND CORRECTION BY CODES 

Decoherence in quantum systems is, according to the current operationalistic point of 
view 



10(1 , induced by coupling the system to its environment. As a result, the states of the 
computer become entangled with the environmental ones, the fact that has a deadly impact 
on quantum superpositions. Formally, decoherence is described by the chain 



Qo 



A) = £?o ® Qo A = U{t)D U\t) — ► ft = Tr He D t , (1) 



where Qo and g t are the computer states (density matrices) at time and t, respectively, 
and £>q is the initial environmental state. Here the leftmost arrow represents enlarging of the 
computer's Hilbert space 7i c by the environmental degrees of freedom, the middle one stands 
for the joint computer-environment evolution (resulting typically in a non-product density 
matrix D t ), and the rightmost arrow reflects our ignorance of the environment, expressed in 
terms of the partial trace over the environmental Hilbert space 7i e . Despite the unitarity of 
the joint evolution D Q — > D t , the full transformation g ~~ > Qt need no longer be unitary. 

During the run of a calculation, the computational state faces a sequence of neatly timed 
unitary transformations. In the meantime, when quantum gates are being readied for the 
next computational step, the quantum registers containing qubits are falling victim to the 
harsh intervention of the environment. Consequently, to study decoherence effects in the 
quantum computer, one should primarily be interested in qubits suspended in the registers 
that are exposed just to the environmental interaction. 

The simplest quantum register contains only one qubit (we denote the single-qubit 
Hilbert space by TC q ). The action of the overall evolution operator U(t) can be formal- 
ized by the following equations 
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U(t)\e t )\0) = \ 9l )\0) + \k)\l), U(t)\ei)\l) = \ Ul )\0) + K>|1>. (2) 

Here |ej) is an orthonormal basis of H e (for simplicity we suppose = J2i w i\ e i)( e i\)i an d 
\gi), \ui) and \rrii) are some unknown environmental states (neither normalization nor 
orthogonality is required) containing t as an implicit variable. The crucial observation [|J, 
leading ultimately to the correcting codes, is that the time-dependence of the evolution 
operator U(t) can be totally embodied in the environmental Hilbert space 7i e . Namely, U(t) 
admits a trivial factorization U{t) = U^{t) ® Q^, where U^it) and (/i = 0, 1, 2, 3) act 
on the environmental and qubit states, respectively. Let = a\0) + The explicit 

form of U^t) and can be deduced from the formula 

U(t)\em = l9i) + 2 lmi) (a|0) + /3|1» + k) ~ K) H0) - /3|1» + 

+ ^^(/3|0) + a|l)) + ^-=^(-/3|0) + a|l)) = 
= 1^)1^) + l&iKIV*) + Iq)^!^) + \di)(-ia y )\^). (3) 

Here a^, a^, and a z are Pauli matrices in the basis |0) and |1). The expression (^) has led 
to a "discrete" classification of quantum errors || into amplitude errors (<J X ), phase errors 
(cr 2 ), and combined amplitude-phase errors {—iu y = u z a x )\ 1 represents no error. In view 
of @, if the three classes of errors can be corrected, then any (even unknown) single qubit 
state can resist perturbations. 

An error decomposition analogous to (|3|) is valid also for n-qubit systems. If U (t) is an 
evolution operator (even non- unitary) on 7i = 7i c ® H c , with 7i c = H® n , then 

^) = E^}W«-<„. (4) 

itn} 

where the sum runs over all n-collections {/ij}™ =1 of indices /ij G {0, 1, 2, 3} and are 
operators on 7i e , corresponding to the respective dynamics in the qubit sector. Here we 
have made use of the shorthand notation 

a k u = 1 ® . . . ® 1 (g) ff M ® 1 g . . . g 1 , (5) 

V V 

1 n— & 
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where a^=(l,a). The decomposition (f|) is unique since it can be inverted to yield 
U{iH}{t) = 2~ n Tr- Wc (?7(t)(j / l ti cr^ 2 . . . ). Now, it is meaningful to say that an error of rank k 
has occurred, if only collections with exactly k non-zero indices \x\ contribute to (|J). 

Error correction by codes relies upon the possibility to encode single logical states |0) and 
|1) into specific superpositions |0) and |1) of compound n-qubit logical states (so-called code- 
words from TC® n ). If the states |0) and |1) are distant enough from each other, where the 
Hamming distance turns out to be an appropriate measure for these purposes, the histories 
of a certain sub-class of errors can be traced back and subsequently remedied. Let k be fixed. 
Then it turns out that for some n satisfying the quantum Hamming and Gilbert- Varshamov 



bounds [11 



£(l><^g(;>', ( 6 , 

there exists a code into n qubits capable of rectifying all errors up to rank k. The above 
inequalities follow from the requirement that the Hilbert space of code-words is large enough 
to allow both, the original information contained in g (cf. formula [T]) as well as the way it 
has been altered by the above sub-class of errors, to be deciphered from the register density 
matrix g t . It will be important for us to observe the asymptotic behaviour of the above 
bounds as n — > oo. Taking at the same time k/n — > x, the formula (P) boils down to 

aTn3 - ln^l - a;) 1 "*] < In 2 < 2sln3 - \n[(2x) 2x (l - 2x) 1 - 2x ]. (7) 

Since both the inequalities are satisfied for a nonvanishing argument (x G [xo,2xo] with 
xq ~ 0.1), we see that the number of errors k that can be controlled grows roughly linearly 
with the size n of the code. 

The repair of defects is standardly thought to be done by means of auxiliary qubits, 
so-called ancillas. Before the correction procedure is initiated, the ancillas are prepared 
in a fixed state \A) (this is important as it implies zero entropy input — disorder can be 
transferred to the ancilla Hilbert space 'Ha). When they are brought into contact with the 
computational qubits, the corresponding type of error is re-printed in the ancilla state and 
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subsequently rectified in the computational sector. Formally, to represent this operation, 



one introduces [^] a recovery operator 1Z, acting on the product space H c ®TLa- As can be 
easily shown, 1Z can be made unitary [ |T2| , |i"3| , which allows for the coherent inclusion of the 
error correction into the computational algorithm. 

III. CORRECTION EFFICIENCY 

It is intuitively clear that the correction procedures function only when the stored com- 
putational state has not departed very far from the initial one. In this section we would 
like to discuss this common supposition quantitatively. In order to do that, let us recall 
|T4| the definition of the fidelity functional ascribed to the channel decoherence + correction 
(represented by the operator 1ZU). Namely, the fidelity functional E#(t) is the probability 
that the initial state passes through the channel intact. If the environment is originally 
in the state Qq, then 

F^t) = Tr[w(t)(P^ ® eg ® P A )U\t)V) P^\. (8) 

Here P^ = is the projector onto the state and, similarly, P A = Without 

being explicitly marked out, the ancilla is not to be affected by the evolution under U(t). 
Consequently, error E(t) of the respective code is defined |12| by E(t) = sup^ E^(t), where 
E^(t) = 1 — F^(t) is the error functional. It is worth noting that as we suppose TZ to be 
unitary, E^ is expressed by the same formula as F^ with only the last projector P^ replaced 
byl-P*. 

Assume now n qubit registers being in the state The registers are exposed to the 
environmental intrusion, described by an interaction hamiltonian V. In the following, we 
shall be concerned with the class of non-contact interactions, i.e., those for which V takes 
the form 

n 3 

v = EE*i.®°i> (9) 

1=1 ,u=l 
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where h are some self-adjoint operators on the environmental Hilbert space (the term with 
\i — is absent in equation (|9]) as it can be incorporated into the free hamiltonian — see 
below). This assumption is justifiable if all inter-qubit communications are mediated by 
some (quasi-) particle external fields. For instance, in the linear ion-trap computer single 
qubits interact by exchanging photons and phonons, so the interaction is indeed a non- 
contact one. 

In the interaction picture, the evolution operator U (t) satisfies the Schrodinger equation 

*jU{t) = V(t)U(t), (10) 

where V(t) is the freely evolved interaction hamiltonian, V(t) = e lH o t y e ~ tH ot ( we se ^ = \y 
As the free hamiltonian is not supposed to induce interaction between qubits (H is the sum 
of environmental and single-qubit terms), the free evolution of ex's in @ can be incorporated 
into the evolution of Tt's. Namely, taking into account that cr l (t) = J2u flf J ,(i) (T l (where f l Vix (t) 
are some functions of time), we can put the operator V(t) to the same form as V, only with 
h], being replaced by £a /^(t)/^*). 

Let the environment be in the state |ej). Then the joint time evolution of |ei)|\&) in the 
interaction picture can be written 

U(t)\ei)\V) = |ei)|tt) -i f ' dsV(s)U(s)\ei)\^) = 

Jo 

= gir(- i /;^v W )'|e < >l*) + 

+ f'd Sl f"ds 2 ... f' k ds M V( 3l )V(, 2 ) . . .V(3 t+1 )U(s t+1 )\e,)\1}, (11) 

Jo Jo Jo 

where T is the time-ordering operator and k is chosen to be precisely the rank of errors we 
intend to correct. The second equality has been obtained by iterating the preceding one k- 
times, while tacitly supposing that U(t) is sufficiently differentiable when applied to |ej)|\l/). 
Now let the ancilla qubits (initially in the state \A)) join the evolution and carry out the 
correction TZ. Since the correction procedure cures all errors up to rank k, the first term 
in ([11]) is thoroughly reverted to the state in the computational sector and, therefore, 
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brings no contribution to the error functional Eq,{t) (due to the presence of 1 — Py). This 
enables us to write 



Ey(t) = ds x ... I ds k+ i / ds[... ds' k+1 
Jo Jo Jo Jo 

Tr[RV( Sl ) . . . V(s k+1 )U(s k+1 )D %A U\s' k+1 )V(s' k+1 ) . . . V{s' x )K\l - P v )], (12) 

with Dy^A = -P* <8> Qq Cg> Pa- We have the first result: // the correcting method is capable 
of remedying all errors up to k, then E^(0) = E^O) = . . . = E% (0) = for all initial 
states 1^). Moreover, the behaviour for short times is given by the relation 

f2k+2 

E*(t) = I ^ EE^W 1 • • • V^Dv.aV 1 '^ . . . V l W{l-P*)) + 0(t 2h+3 ), (13) 

where the sums run over all ordered k + 1-subsets of {1, 2, . . . , n}, and V 1 = ELi h l (0) <8> cr l . 
The formula fll3f) is obtained by substituting V(s) ~ V"(0) and U(s) ~ 1 for small times 
s in (p!2[ ). As the errors incurred in (jT3"|) by the chains of V's are of rank k + 1, they are 
not all likely to be corrected for a generic environmental interaction, unless the correction 
method encompasses also higher order errors. Thus, E^ h+2 \o) ^ and (|T3|) describes the 
true behaviour of Ey for short times. 

We see that, indeed, the more enhanced codes are applied the slower the errors escalate. 
It is worth noting that the formula (|T3|) generalizes the standard treatment of the quantum 
watch-dog effect. In the latter, no error correction is employed, therefore k — 0, which is 
in accord with the obtained time dependence (E^(t) ~ t 2 ). However, as opposed to the 
watch-dog effect, quantum error correction by codes requires no knowledge of the initial 
state. 



IV. STABILIZATION BY CORRECTION 

We have seen in the previous section that the ancilla based correction brings about 
polynomial slow-down of error propagation, provided the mild condition (|9|) upon the form 
of the interaction is satisfied. In practice, however, this does not imply that decoherence 
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can be stopped from penetrating into the computer. Namely, as already mentioned, since 
the frequency at which the rectification is repeated cannot be made arbitrarily large, even 
the encoded quantum information in the register subjected to periodic corrections decays 
exponentially with time, with only the decay rate reduced. In view of these remarks, the 
following concept of stabilization is more appropriate: A correction procedure, depending on 
a discrete parameter n, qualifies to be stabilizing, if there is a range of times for which the 
probability of errors can be made arbitrarily small by varying the parameter n. The latter 
parameter represents the comprehensiveness of the correction method and is exemplified by 
the length of code-words in our case. 

We would like to clarify when a computer register can be stabilized in the above sense. 
In particular, we prove the following statement: If the decohering interaction has the form 
(^), with the operators h 1 ^ uniformly bounded, then there exists an error procedure stabilizing 
quantum computer against error proliferation. It should be noted that we make no particular 
assumptions on the free hamiltonian, apart from the natural self-adjointness requirement 
ensuring the existence of e~ lHot . For the proof, notice that in this setting all the above 
formulas can be given a good mathematical meaning. Consequently, we can derive the 
inequality 



by directly utilizing the boundedness assumption in (|l^) (|| ■ || is the standard operator 
norm). If \\h l \\ < C for all \i and I, then, in the regime k,n —>■ oo, we get 



where Xq stands for the lower bound on the asymptotic value of k/n (see (|7|)). Hence, if 
t < Xo/Ce, the error E(t) tends to zero exponentially fast as n increases. Consequently, 
by choosing a sufficiently enhanced encoding, the fidelity of the decoherence+correction 
channel outstrips every bound 1 — 5 (with 5 > 0) and the quantum information is, indeed, 
well stabilized. 




(14) 



E(t) < 0(1) tCe 




(15) 
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Unfortunately, it turns out that the genuine environmental interaction cannot be rep- 
resented by a bounded operator. Nevertheless, the high-energy environmental states are 
usually strongly suppressed due to low temperature T. Since the characteristic scale of the 
interaction energy h 1 ^ is of order k&T, one is led to the rough estimate t < tx ~ l/feT 
for the regime of applicability of the quantum error correction. Hence, the threshold of the 
thermal regime tx (cf. ||) probably sets an absolute bound beyond which no correction is 
of substantial help. 

The rigorous treatment of the unbounded case in this generality falls beyond the current 
level of mathematical quantum theory. It is not so difficult to control the formula (0) if 
certain plausible assumptions on the interaction H + V can be made. In particular, one 
goes from the description in terms of the Dyson series to the description by the Born series 



| 15| . It is not clear, however, whether it is possible to derive an analogue of the bound 
(f[4j)) in particular, with the factor (k + l)! 2 in the denominator (notice that this factor is 
absolutely necessary to compensate the natural extensive behaviour of the operator V). We 
plan to address these topics in a future work. 

The restriction to non-contact interactions is, actually, not essential. Our treatment 
can be applied also to contact hamiltonians (with multiple-a terms), leading, however, to 
correspondingly weaker results concerning the correction efficiency. 
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